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ABSTRACT 


Rotationally  mediated  selective  adsorption  scattering  resonances 
are  used  to  make  an  experimental  and  theoretical  study  of  the  laterally 
averaged  interaction  potential  between  HO  and  a  weakly  corrugated  system, 
Ag(lll).  The  experimentally  observed  resonances  determine  the  vibrational 
levels  of  the  HD/ Ag( 111)  physisorption  potential  as  a  function  of  bound 
rotational  state.  These  vibrational  levels  show  J-dependent  shifts  due 
to  the  orientational  anisotropy  of  the  potential .  Exact  quantum  scattering 
calculations  using  a  full  laterally  averaged  potential  of  the  form 

VQ{ 2,0)  =  vQ(*z)  Cl  +  6P2  (cose) 3 

have  been  carried  out  to  obtain  rotationally  inelastic  transition  probabil ities . 
Experimental  and  theoretical  resonance  energies  are  compared  for  two  forms  of 
v0(z),  a  Morse  and  a  variable  exponent  potential,  as  a  function  of  6,  and  are 
found  to  be  very  close  to  the  first  order  perturbed  energies  of  a  free  rotor 
in  bound  states  of  vQ(z),  Both  potential  forms  give  equally  good  fits  to  the 
data,  yielding  an  optimum  value  of  the  asymmetry  parameter, 6 s -0.05.  The 
determination  of  8  is  relatively  insensitive  to  small  changes  in  the  vQ(z) 
well  depth. 
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1.  INTRODUCTION 

The  scattering  of  atomic  and  molecular  beams  from  well -characterized 

single  crystal  surfaces  is  attracting  considerable  experimental  and 
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theoretical  interest  at  this  time  ’  .  One  of  the  primary  objectives  of 
these  studies  is  to  determine  the  nature  of  the  interaction  potentials 
which  exist  between  a  wide  variety  of  incident  atomic  and  molecular 
particles  and  different  classes  of  clean  or  adsorbate  covered  crystalline 
surfaces  including,  for  example,  metals,  alkali  halides,  semiconductors, 
graphite  and  rare  gases.  Knowledge  of  these  gas-surface  potentials  is 
central  to  our  understanding  of  many  important  heterogeneous  processes 
since  they  determine  the  clynamlcs  and  energetics,  i.e-,the  time  evolution 
of  these  phenomena  on  the  molecular  level.  In  addition,  very  accurate 
Information  on  surface  electron  charge  density  distributions  and  surface 
structures  can  also  be  extracted  from  diffraction  experiments,  provided 
that  the  electrostatic  terms  contributing  to  the  overall  gas-surface 
potential  are  known.  For  helium-surface  Interactions  a  theoretical 
framework  which  analytically  relates  surface  electron  charge  densities  to 
the  helium-surface  potential  is  now  in  hand"3  ,  and  has  been  used  in  con¬ 
junction  with  first  principles  surface  charge  density  calculations  to  obtain 
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structural  information  from  diffraction  data  . 

Recent  attention  in  this  field  has  been  focused  on  the  study  of 

Interaction  potentials  between  helium,  atomic  hydrogen,  molecular  hydrogen, 

1  2 

and  several  highly  corrugated  clean  surfaces,  ’  and  physisorbed  layers  of 
6  7  8 

xenon  on  graphite  ’  and  silver  .  One  of  the  primary  reasons  for  studying 
these  particular  systems  is  the  large  non-specular  diffraction  peaks  which 
they  show.  Consequently,  very  strong  diffractive  selective  adsorption 
scattering  resonances  can  be  observed  for  these  systems.  Such  diffractive 
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selective  adsorption  resonances,  which  were  first  observed  in  1930  by 
9 

Estermann  and  Stern  and  subsequently  explained  by  Lennard-Jones  and 
Devonshire,1*1  are  extremely  useful  as  they  can  be  used  to  determine  the 
eigenvalues  of  the  laterally  averaged  particle-surface  interaction 
potential.  The  resulting  eigenvalue  spectrum  can  then  be  used  to  determine 
the  shape  (but  not  absolute  location)  of  the  laterally  averaged  potential 
with  respect  to  the  surface  plane.by  either  fitting  the  spectrum  to  that 
of  a  model  potential,  or  by  a  procedure  analogous  to  a  Rydberg-Klein-Rees 
(RKR)  potential  inversion  .  In  some  special  instances  splittings  are 
observed  when  two  diffractive  selective  adsorption  resonances  are  mixed 
by  the  periodic  part  of  the  gas-surface  potential.  When  this  occurs 
Information  on  the  higher  Fourier  coefficients  of  the  periodic  potential 
expansion  can  be  also  be  obtained  from  analysis  of  the  band  structure  of 
the  resonances15. 

In  our  laboratory  we  are  currently  attempting  to  extend  the  above 
studies  to  include  the  detailed  aspects  of  molecular  interactions  with 
relatively  smooth  metallic  surfaces.  Several  scattering  experiments  and 
complementary  close-coupled  scattering  calculations  have  now  been  carried 
out  in  order  to  quantitatively  determine  both  the  isotropic  and  anisotropic 
(orientation  dependent)  components  of  the  laterally  averaged  molecular 
hydrogen  -  Ag( 111)  potential.  The  experimental  studies  include  determination 
of  the  bound  physisorption  levels  of  HD  on  silver  as  a  function  of  rotor 
state, using  our  newly  developed  technique  of  rotational ly  mediated  selective 
adsorption  (hereafter,  RMSA) 15 ;  determination  of  the  physisorption  levels 
for  and  D£  on  silver  using  diffractive  selective  adsorption  (hereafter, 
DSA)  scattering  resonances;  and  determination  of  the  rotationally  inelastic 
scattering  cross-sections  as  a  function  of  collision  energy  for  Hg,  and 


HD  reflecting  from  silver. 


The  motivation  for  making  the  above  measurements  is  quite  straightforward. 

The  isotropic  "zeroth -order"  physisorption  potential  can  best  be  determined  by 
mapping  out  the  J=0  bound  state  resonances  for  H2  and  D2  using  DSA.  Once  this 
is  determined, DSA  measurements  for  J=1  H2  and  Dg,  and  RMSA  measurements  for 
HD( J>0)  can  be  used  to  obtain  the  anisotropic  component  of  the  molecule-surface 
potential.  The  basis  for  this  is  that  the  anisotropy  of  the  laterally  averaged 
potential  gives  rise  to  J  and  m  dependent  shifts  of  the  observed  bound  state  eigen¬ 
values.  For  both  DSA  and  RMSA  measurements  the  eigenvalue  spectrum  is  given  by  first 
order  perturbation  theory^.  Analysis  shows  that  DSA  measurements  are  dependent 
upon  both  the  J  and  m  states  of  the  bound  rotor,  whereas  RMSA  measurements  on  a 
weakly  corrugated  surface  show  J-dependence  only,  since  m  is  conserved  for  a 
flat  surface.  (The  incident  HD  ensemble  used  in  these  studies  was  predominantly 
in  the  |J,m>  =  |0,0>  state).  An  estimate  of  the  anisotropic  component  can  then 
be  arrived  at  by  comparing  the  experimental  resonance  energies  with  those  obtained 
from  the  close-coupled  calculations.  The  rotationally  inelastic  H2  and  D2  cross- 
sections  are  also  quite  sensitive  to  the  anisotropy  of  the  interaction  and  will 
utimately  be  used  as  an  independent  test  of  the  potential  derived  from  the  RMSA 
and  DSA  studies . 

In  this  paper  we  present  one  aspect  of  this  comprehensive  study,  namely 
the  experimental  and  theoretical  results  for  the  RMSA  studies  of  HD  inter¬ 
acting  with  Ag(lll)  as  a  function  of  bound  rotational  state.  A  brief 
review  of  the  RMSA  process  is  presented  in  the  results  section  of  the  paper. 

This  work  follows  our  earlier  study  of  HD/Pt(  111)  interactions^  and  our 

2 

subsequent  demonstration  that  Wigner  (L  )  R-matrix  scattering  techniques 

can  be  used  to  investigate  the  nature  of  molecule-surface  bound  states^. 

18 

Similar  quantum  calculations  have  also  been  recently  reported  by  Schinke 
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who  treated  the  rotational  degree  of  freedom  exactly  within  the  diffraction 
sudden  approximation  .  We  feel  that  Ag(lll)  is  a  nearly  ideal  surface 
for  these  studies  due  to  its  Inert  chemical  properties  with  respect  to 
molecular  hydrogen.  This  has  allowed  us  to  cryogenical ly  cool  the  substrate 
in  order  to  diminish  the  diffuse  Inelastic  component  and  Deby e-Waller 
attenuation  of  the  scattered  molecules,  relative  to  that  which  was  achieved 
with  Pt( 111) . 

It  should  be  emphasized  that  the  RMSA  results  presented  in  this  paper 
have  been  used  to  determine  the  anisotropic  component  of  the  hydrogen/Ag(lll) 
potential  by  building  upon  the  "zeroth-order"  isotropic  potential  which  was 
recently  obtained  in  our  DSA  studies  of  Hg  and  on  Ag(lll)  .  It  is  our 
belief  that  the  unique  and  complementary  information  obtained  by  combining 
RMSA  and  DSA  measurements  will  enable  the  determination  of  hydrogen/surface 
interaction  potentials  for  a  wide  variety  of  substrate  materials. 

II.  APPARATUS  AND  EXPERIMENTAL  PROCEDURE 

The  scattering  apparatus  used  in  these  studies  can  be  schematically 

portrayed  as  consisting  of  four  main  sections  which  are  the  beam  source 

and  vacuum  manifold,  detector,  crystal  assembly  and  associated  surface 

analysis  equipment,  and  data  acquisition  electronics.  Each  of  these 

sections  will  be  briefly  reviewed  before  describing  our  experimental 

procedures.  A  detailed  description  of  the  apparatus  has  been  presented 
21 

previously  . 

The  HD  needed  for  these  experiments  was  generated  catalytically  by 
passing  an  8:1  D^  mixture  through  a  600  K  stainless  steel  column  which 
was  packed  with  magnesium  granules.  The  HD  generator  is  operated  on-line, 
and  has  worked  reliably  for  well  over  one  year.  The  H^:  0 ^  ratio  was  chosen 


to  be  8:1  in  order  to  make  the  best  compromise  between  producing  a  beam  of 

relatively  low  velocity  dispersion  and  high  HD  intensity.  The  isotopically 

equilibrated  effluent  of  this  column  (nominally  4:1  is  passed 

through  an  in-line  liquid  nitrogen  cooled  molecular  sieve  trap  which  serves 

to  remove  any  moisture  coming  from  the  HD  generator  or  from  the  feed  gas  itself. 

This  trap  was  necessary  as  most  of  the  experimental  data  was  collected  with 

the  crystal  held  at  cryogenic  temperatures.  The  beam  source  contained  a 

20  micron  pinhole  and  could  be  ohmically  heated  or  cryogenically  cooled  with 

liquid  nitrogen.  Three  different  beam  energies  were  used:  Eb=40.6  meV 

(PN  =  150  PSI,  Av/v  =  5.0%  FWHM) ,  Eb=113.6  meV  (PN  =  250  PSI,  Av/v  =  5.2%) 

and  Eb=153.0  meV  (P^  =  300  PSI,  Av/v  =  6.0%).  An  el ect reformed  skimmer 

22 

was  used  to  ensure  a  high  terminal  Mach  number  for  the  HD  expansion  . 

The  nozzle  beam  is  differentially  pumped  in  three  stages  in  a  non-bakeable 
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source  chamber,  with  typical  operating  pressure  being  5X10  ,  1X10  ,  and 

5xl0~7  torr  for  the  regions  before  the  skimmer,  before  the  beam  collimator, 

and  containing  the  chopper,  respectively.  The  main  UHV  scattering  chamber 

is  separated  from  the  three  source  chamber  regions  by  a  bakeable  buffer  region 

_g 

whose  pressure  rises  to  5x10  torr  when  the  beam  is  on.  This  buffer  section 
is  separated  from  the  chopper  section  by  a  sliding  gate  valve  which  contains 
the  0.5mm  diam.  beam  defining  aperture.  Supported  by  a  400  1/sec  ion 

pump  and  a  titanium  sublimation  pump,  the  crystal  chamber  has  a  base  pressure 
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of  ^  7x10  torr,  which  rises  as  high  as  1x10  torr  when  the  beam  is  on. 

The  incident  beam  is  collimated  to  an  angular  divergence  of  0.1°,  which 

corresponds  to  a  0.5mm  diam.  beam  spot  at  the  crystal. 

Particle  detection  is  accomplished  with  a  doubly  differentially  pumped 

quadruple  mass  spectrometer,  which  rotates  about  a  horizontal  axis  in  the 
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plane  containing  the  incident  beam  and  the  surface  normal.  The 
electron  bombardment  ionizer  is  located  14.45  cm  from  the  sample,  and 
is  collimated  to  view  the  crystal  with  1.2°  resolution.  The  detector 
can  be  rotated  manually  or  under  computer  control. 

The  Ag(lll)  crystal  was  cut  from  a  99.999%  boule  and  polished  at 
the  Cornell  Materials  Preparation  Laboratory.  It  was  subsequently  re¬ 
polished  in  our  laboratory  with  very  fine  cerium  oxide.  This  was  followed 
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by  electro-chemical  etching  with  a  saturated  chromic  acid- 10%  HC1  solution  . 

The  residue  from  this  treatment  was  removed  with  a  10%  HNO^  solution.  The 
orientation  of  the  crystal  was  checked  by  Laue  X-ray  back  reflection  and 
was  found  to  be  within  of  the  (111)  plane.  The  underside  of  the  crystal 
disc  was  spotwelded  to  two  1mm  thick  Pt  wires,  which  were  attached  through 
two  tungsten  rods  to  the  manipulator.  A  Chromel -Alumel  thermocouple  was  also 
spotwelded  to  the  back  of  the  crystal.  The  crystal  can  be  resistively  heated 
or  conductively  cooled  through  the  Pt  support  wires.  Temperature  regulation 
is  achieved  with  a  home  built  circuit  which  controls  the  current  source. 

The  crystal  can  be  translated  in  three  mutually  perpendicular  directions, 
tilted  out  of  plane,  and  rotated  about  two  orthogonal  axes  to  set  its 
azimuthal  and  polar  orientation  with  respect  to  the  incident  beam. 

After  bakeout,  the  crystal  was  subjected  to  several  cycles  of  Ar  ion 
sputtering  (2  uamps  at  500  V  for  10  minutes)  until  no  appreciable  contamination 
could  be  detected  with  a  single-pass  CMA  Auger  spectrometer.  The  crystal 
was  then  annealed  at  600  K  for  one  hour  to  remove  sputter  damage.  The 
success  of  this  annealing  procedure  was  verified  by  the  fairly  uniform  He  and 
reflectivities  which  were  measured  along  several  points  of  the  crystal. 
Debye-Waller  corrected  reflectivities  ( i  .e.,  extrapolated  to  T<.  =  0  K)  of 
85%  were  observed  for  incident  at  tK  =  70°,  50°,  and  30°  with  respect  to 


the  surface  normal. 

Data  acquisition  and  instrument  control  is  carried  out  with  a  PDPll 
minicomputer  which  is  interfaced  to  the  scattering  apparatus  via  CAMAC. 
Time-of- flight  measurements  are  made  with  a  255  channel  multiscaler 
which  has  a  minimum  dwell  time  of  0.25  usee/ channel .  Molecular  beam  gating 
is  carried  out  with  a  15.24  cm  diam.  chopper  which  contains  two  narrow 
slots  for  velocity  analysis  (0.5  mm  width)  and  two  50%  duty  cycle  patterns 
for  angular  distribution  measurements.  Angular  distributions  are  obtained 
with  two  32  bit  scalers  which  are  synchronized  to  the  chopper  in  order  to 
collect  signal  plus  background  (slot  open)  and  background  only  (slot  closed) 
information.  RMSA  signals  were  collected  by  a  Keithley  electrometer, 
demodulated  by  a  lock-in  amplifier,  passed  to  a  digital  voltmeter,  and  finally 
read  into  the  computer. 

In  the  next  section  two  types  of  experimental  data  are  presented:  an 
HD  angular  distribution  and  several  RMSA  scans.  Each  experiment  was  always 
preceeded  and  immediately  followed  by  a  beam  characterizing  time-of- fl ight 
measurement.  The  angular  distribution  was  obtained  by  scanning  the  detector 
in  0.2°  increments  for  fixed  orientation  of  the  crystal  with  respect  to  the 
incident  be-m.  The  RMSA  scans  were  taken  by  measuring  the  maximum  intensity 
of  each  outgoing  rotationally  elastic  or  inelastic  diffraction  peak  as  a 
function  of  incident  angle.  In  practice,  data  was  collected  for  only  one 
diffraction  peak  at  a  time  as  0.  was  varied.  This  allowed  us  to  alternately 
move  the  crystal  polar  angle  by  ca.  0.2°,  and  then  scan  the  detector  through 
the  same  peak  with  the  computer  determining  the  angle  of  peak  reflected 
intensity  by  recording  in  real  time  the  signal  intensity  versus  detector  angle 
The  peak  heights  determined  in  this  manner  include  contributions  from  both 
the  coherently  scattered  particles  and  the  diffuse  inelastic  background. 


8 


III.  EXPERIMENTAL  RESULTS 

The  RMSA  scans  are  based  upon  the  reflected  particle  angular 
distributions  for  a  large  range  of  incident  angles.  It  is  therefore 
appropriate  to  begin  this  section  by  presenting  a  typical  HD  angular 
distribution  for  one  (arbitrarily  chosen)  incident  angle.  This  should 
help  clarify  the  kinematics  which  are  an  integral  part  of  the  RMSA  data. 

Such  an  angular  distribution  is  shown  in  Figure  1. 

The  discrete  peaks  seen  in  this  figure  are  due  to  elastic  and 
rotationally  inelastic  diffraction  from  the  relatively  flat  and  well-ordered 
Ag(lll)  surface.  Similar  scattering  patterns  have  been  previously  reported 
for  hydrogen  isotopes  reflecting  from  MgO(OOl)2^,  LiF(OOl)2^,  and  Pt  ( 111)  16»^6.27 
The  initial  ensemble  of  HD  molecules  is  predominantly  in  the  |J,m>  a  |0,0> 
state.  The  low  population  of  higher  rotational  states  is  due  to  both  the 
high  quality  of  our  supersonic  jet  expansion,  and  to  the  absence  of  the  aJ=2 
selection  rule  for  HD.  Rotationally  inelastic  peaks  are  labelled  where 

these  labels  represent  the  initial  and  final  rotational  states  of  the  reflected 
particles.  These  inelastic  transitions  are  associated  with  the  (00), (10)  and 
(10)  diffraction  peaks.  Their  final  scattering  angles  are  given  by  the 
momentum  and  energy  conservation  constraints  which  govern  the  scattering 
kinematics.  These  are  (i)  conservation  of  parallel  momentum 


■  «,  ♦ «  <»> 

and  (ii)  conservation  of  total  energy 

=  ,£f2  +  (2M/ii2)  CE(Jf)  -  E(Ji)3  (2) 


iW  iiKrtd 


where  J^.  and  are  the  initial  and  final  wave  vectors,  M  is  the  mass  of  HD, 
E(Jf)  -  is  the  change  in  rotational  energy,  J<i  and  ^  are  the 

components  of  ^  and  J<^  parallel  to  the  surface,  and  ^  is  a  reciprocal  lattice 
vector.  The  above  equations  do  not  include  phonon  interactions  which  are 
weak  processes  for  the  collision  energies  and  surface  temperatures  used  in 
this  study.  The  peak  assignments  shown  in  Figure  1  have  been  independently 
verified  with  time-of-flight  velocity  measurements. 

The  lower  portion  of  Figure  1  is  an  Ewald  representation  of  the  scattering 
kinematics.  This  provides  us  with  a  very  useful  visualization  of  the  scatter¬ 
ing  processes  which  are  occurring.  The  radii  of  the  five  semicircles  represent 
the  magnitudes  of  the  final  k-vectors  associated  with  the  different  transitions. 
Energy  conservation  requires  that  the  final  elastic  and  rotationally  inelastic 
peaks  fall  on  the  appropriate  semicircular  arc,  as  given  by  Eq.  (2).  In 
addition,  parallel  momentum  conservation,  Eq.  (1),  requires  that  each  peak  fall 
on  either  the  specular  or  £  f  0  diffraction  rods  (Brillouin  zone  boundaries). 

The  intersection  of  the  semicircles  and  vertical  rods,  therefore,  provide  a 
graphical  determination  of  the  allowed  final  scattering  angles.  The  indicated 
angles  are  those  appropriate  for  a  monoenergetic  HD  beam  having  |k^|  =  12.7 
and  are  in  excellent  agreement  with  the  experimental  results,  justifying 
neglect  of  phonon  interactions  in  the  above  discussion.  The  angular 
width  and  exact  shape  for  each  diffraction  peak  can  be  determined  by 
carrying  out  a  full  convolution  of  the  true  data  with  the  incident  beam's 
velocity  distribution  and  other  instrumental  resolution  parameters.  The 
reflected  peaks  have  angular  profiles  which  are  not  significantly  broadened 
by  surface  domain  size  effects,  indicating  that  the  Ag(lll)  surface  is  well 
ordered.  Scattering  experiments  have  also  been  carried  out  along  the  < 1 12> 
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and  <101>  directions:  no  apparent  change  is  seen  in  the  inelastic  transition 
probabilities  associated  with  the  specular  channel. 

Two  striking  results  are  evident  from  the  angular  distribution  shown 
in  Figure  1.  These  are  the  extremely  high  rotationally  Inelastic  cross- 
sections  for  HD  relative  to  H ^  and  Dg,  and  the  presence  of  only  weak  diffraction 
for  G  t  0.  The  first  observation  can  be  attributed  to  the  offset  of  the  HD 
center-of-mass  from  its  geometric  center.  The  second  observation  indicates 
that  the  HD/Ag( 111)  potential  has  low  corrugation.  This  will  allow  us  to 
treat  Ag(lll)  to  first  approximation  as  a  flat  surface  in  the  scattering 
calculations  presented  later. 

In  Figure  2,  RMSA  scans  (again,  rotationally  mediated  selective  adsorption) 

are  shown  for  HD  incident  along  the  <10 1>  and  <112>  azimuthal  directions.  Each 

of  these  scans  consist  of  three  Independent  data  collection  runs  in  which  the 

intensities  (peak  heights)  of  the  J=0-*0,  0->-l,  and  0-*-2  transitions  have  been 

plotted  verses  E  ,  the  perpendicular  component  of  the  incident  beam  energy. 
zi 

The  actual  data  were  collected  by  varying  the  incident  polar  angle,  in  ca. 

0.2°  increments  in  the  manner  described  earlier  under  experimental  procedures. 

The  sharp  undulations  seen  In  Figure  2  are  due  to  the  phenomena  of 
16  17  18 

RMSA  ’  ’  .  To  a  first  approximation,  these  bound  state  scattering 

resonances  arise  when  W,  /2M,  as  calculated  from  Eqs .  (1)  and  (2)  for 

zf 

a  given  translation-rotation  energy  exchange  process,  is  negative  and 
precisely  matches  a  bound  energy  level  of  the  HD/Ag( 111)  physisorption 
potential.  Resonance  dips  are  observed  in  a  J=0-+n  transition  when  the  J=0+n+l 
transition  "deposits"  the  molecule  into  a  bound  level  of  the  gas-surface  potential 


* 
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(i.e.,  J=0-*-n  +  1  leads  to  a  closed,  rotational  scattering  channel). 

(Hereafter,  we  shall  use  to  denote  when  Is  the  rotational 

label  for  a  bound  molecule).  As  seen  in  Figure  2,  the  J^n  +  1 

resonances  most  strongly  perturb  the  J=0-+n  transitions,  causing  pronounced 

minima  to  appear.  Smaller  complementary  peaks  can  also  be  observed  In 

J=0-*-n  -  1  transitions.  It  therefore  appears  that  the  Jfa=n  +  1  resonance 

states  decay  primarily  to  states  other  than  the  nearest  open  channel,  J=0-*-n. 

Figures  2  and  3  have  intentionally  been  plotted  with  Ez  rather  than 

along  the  horizontal  axis  in  order  to  demonstrate  that  the  resonance  positions 

are  independent  of  both  azimuthal  crystal  orientation  and  incident  momentum 

parallel  to  the  crystal  surface.  This  requirement  of  the  proposed  RMSA 

process  is  clearly  met,  since  in  Figure  2  the  dip  positions  are  invariant 

with  respect  to  azimuthal  orientation,  and  in  Figure  3  the  dip  positions 

are  shown  to  be  invariant  for  three  different  incident  beam  energies  (dip 

positions  are  unchanged  even  though  equivalent  E  values  occur  at  different 

zi 

ei  settings).  In  Figure  3  it  can  also  be  seen  that  the  effect  of  increasing 

the  total  collision  energy  with  E  held  constant  (i.e.,  increasing  the 

zi 

parallel  component  of  incident  momentum)  is  to  increase  the  likelihood  of 
rotational  excitation  relative  to  elastic  scattering. 

In  Table  1  we  present  our  main  experimental  results,  which  are  the 
resonance  energies  for  J^=l , 2  and  3  as  determined  from  the  positions  of 
minima  in  the  J=0-»0,  0-+1,  and  0-+2  RMSA  scans,  respectively.  These 
resonance  energies  were  arrived  at  by  averaging  the  data  taken  with 
different  beam  energies  along  both  the  <1I2>  and  < 10 1>  directions.  Appearing 
next  to  these  values  are  the  approximate  bound  vibrational  levels  of  the 
HD/Ag( 111)  physisorption  potential  as  a  function  bound  rotor  state.  These 


energy  levels  were  calculated  by  subtracting  free  rotor  HD  energies 
from  the  appropriate  resonance  energies. 
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The  discrepancies  between  the  results  for  different  states  imply 

that  HD  behaves  as  a  slighly  hindered  rotor  in  the  bound  levels.  It  Is 

apparent  that  there  is  a  small  but  consistent  shift  in  the  energy  of  a 

given  bound  level  as  a  function  of  J^,  a  shift  which  is  well  outside 

the  experimental  uncertainty  of  our  resonance  energies.  The  small  size 

of  these  shifts  indicates  that  the  hindering  is  relatively  weak.  This  is 

independently  confirmed  by  recent  low  energy  electron  loss  experiments 

29  30 

involving  hydrogen  physisorbed  on  Ag(lll)  and  Cu(OOl)  . 

The  origin  of  these  shifts  and  the  way  in  which  we  can  use  them  to 
extract  the  anisotropic  part  of  the  laterally  averaged  HD/Ag( 111)  potential 
will  be  described  in  the  following  section. 


IV.  THEORETICAL  BACKGROUND 


If  inelastic  processes  between  the  surface  and  incident  molecules  are 
neglected,  it  is  straightforward  to  write  down  the  coupled  equations  for 
scattering  of  a  diatomic  molecule  from  a  corrugated  surface. 

Assuming  periodic  corrugation  and  treating  the  molecule  as  a  rigid 
rotor,  the  Hamiltonian  for  motion  of  the  molecule  is 


H  =  =ft'iki'ies,neae 


£V>^)eL*;S 


Here  z  is  the  normal  distance  from  the  surface  to  the  molecule  center  of  mass, 

Jl  °  (x.y)  is  a  surface  lattice  vector,  and  £  a  reciprocal  lattice  vector.  0  and 
p  are  the  polar  and  azimuthal  angles  of  the  molecular  axis  relative  to  the  z 
axis.  The  molecular  wave  function  may  then  be  expanded  in  products  of 
spherical  harmonics  and  plane  waves  of  the  two  dimensional  surface  net: 


(4) 


't'fee.MV  £  X  W  V <•.*)*' 


c6W$)*& 


where  j<  is  the  incident  parallel  component  of  the  wave  vector.  This  leads 


"here  =  ^JifE ~(K  +  (rf-  T(7+|)^-  <£) 


Since  diffraction  is  very  much  weaker  than  rotational  excitation  for 
the  asymmetric  rotor-metal  system,  it  is  valid  to  neglect  coupling  terms 
for  which  £  -  £7  0  In  studying  the  rotatlonaily  inelastic  spectra.  This 
assumption  of  a  flat  surface  is  a  significant  simplifying  feature,  since  the 
number  of  open  channels  for  HD  at  beam  energies  of  100  meV  is  reduced  from 
about  100  to  4.  The  potential  VQ(z,Q,<p)  is  then  assumed  to  be  independent 
of  <j>,  ensuring  conservation  of  the  azimuthal  quantum  number  m. 

The  potential  to  be  used  for  interpreting  the  rotationally  inelastic 
spectra  is  then  Vo(z,0).  For  a  symmetric  rotor,  an  expansion  of  this  in 
Legendre  polynomials  of  cose  reduces  to  a  sum  over  even  terms  only. 

Y<*»»  £  v*)P(c«e> 

4*o, a,*., .  x 

Here  z’  is  the  normal  distance  from  the  surface  to  the  center  of  mass,  which 
is  the  geometric  center  of  the  molecule  in  this  case. 
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In  the  previous  quantum  calculations  of  scattering  of  diatomic  molecules 

from  solid  surfaces,  the  expansion  is  truncated  after  the  second  term,  by 

17  19  31 

analogy  with  atom-molecule  potential  surfaces  *  *  .  This  Is  exact  for 

asymptotically  large  z,  where  the  molecule-surface  interaction  is  dominated 

32 

by  the  attractive  polarization  force  .  At  small  z,  however,  the  repulsive 

potential  arising  from  overlap  of  electronic  densities  of  the  molecule  and 

of  the  surface  is  dominant.  It  has  not  yet  been  established  whether  higher 

terms  are  present  here  or  whether,  if  present,  it  is  warranted  to  neglect 

them.  However,  since  we  aim  to  provide  the  simplest  possible  parameterization 

of  the  potential  VQ(z, 8)  which  adequately  represents  the  rotationally  inelastic 

selective  adsorption  features,  we  regard  it  sufficient  at  this  stage  to  retain 

31 

only  1=0,2  In  Eq.  (7).  We  further  follow  the  form  of  previous  calculations 
by  making  v’2(z)  =  6vq(z),  where  vQ(z)  is  the  laterally  averaged  potential 
for  interaction  of  a  structureless  particle  with  the  surface.  The  potential 
is  then  entirely  determined  by  8  and  vQ(z). 

The  corresponding  potential  for  a  chemically  equivalent  (isotopic) 
asymmetric  rotor  is  obtained  by  transforming  VQ(z',0)  from  the  geometric 
center  of  the  diatom  to  the  displaced  center  of  mass; 

V  « ©)  =  V0  (2)  [  I  +  p  ?2  (  cos  ©)]  <8a> 

=  v0  (-a)  +  v'  Cz) P  (cose)  (8b) 

*  a 

Here  z1  is  the  normal  distance  from  the  surface  to  the  geometric  center  of  HD, 
and  z  is  the  normal  distance  from  the  surface  to  the  center  of  mass: 

z'  =  z  +  <5  cose  (9) 

where  6  is  the  displacement  of  the  center  of  mass  from  the  geometric  center. 

For  HD,  6=0.233ao . 
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The  laterally  averaged  potential  vQ(z)  is  separated  out  to  obtain  a 
decoupled  zeroth-order  Hamiltonian, 


The  full  Hamiltonian  is 


with 


(11) 

(12) 


H  -  H  +  V 

•V  *  «» 

Eigenfunctions  of  H°  will  be  referred  to  as  a  bound- vibration- rigid  rotor 
states  (BVRR  states),  following  ref.  (17).  They  are  direct  products 
of  eigenfunctions  of  vQ{z),  which  are  exactly  analogous  to  the  bound  states 
probed  in  atomic  diffractive  selective  adsorption,  and  rigid  rotor  functions. 

If  the  resonant  states  were  unhindered  rigid  rotors  in  bound  states  of 
vQ(z),  the  resonance  energies  would  be  equal  to  the  BVRR  energies.  However, 
the  BVRR  states  are  not  eigenvalues  of  the  full  Hamiltonian:  the  couplingN/" 
Eq.  (11),  will  shift ‘the  resonance  energies  from  the  BVRR  energies.  The 
extent  of  the  shift  is  an  indication  of  the  amount  of  interaction  between 
rotational  and  vibrational  motions  in  the  quasibound  state.  It  was  recently 
shown1^  that  for  a  potential  of  the  form  In  Eq.  (8)  appropriate  to  HD  on 
metals,  resonances  in  the  inelastic  transition  probabilities  occur  at  energies 
very  close  to  the  first  order  perturbed  BVRR  energies. 

Cl  2  -0  (13) 

“  6n  T(T+0 

c{JL 


Here  Ec  is  the  calculated  resonance  energy,  E^  j  is  the  BVRR  energy  plus  the 
first  order  correction  from  the  potentiaTN/^  Eq.  (12),  e  is  an  eigenvalue 
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of  vQ(z)  and  |nJ>  is  a  BVRR  eigenfunction, with  n  denoting  the  vibrational 
state,  and  J  the  rotor  state. 

We  used  the  same  method  of  solving  the  close  coupled  equations  here 

as  in  ref.  (17),  namely  the  single  R-matrix,  generating  BVRR  states,  energies 

E°  ,  and  £'  ,,  and  rotationally  inelastic  transition  probabilities.  We  did 

n  J  n  J 

this  for  two  different  forms  of  potential  vQ( 2) ,  with  a  range  of  values  of  6  in  each 

case.  The  R-matrix  calculation  employs  a  basis  of  the  direct  product  of  free 

rotor  states  with  a  translational  basis  in  which  the  bound  vibrational  states 

are  expanded  on  a  finite  Inner  region.  Sine  functions  were  used  for  the 

translational  basis.  Increasing  the  size  of  this  region  and  the  size  of  the 

translational  basis  allows  the  bound  vibrational  states  to  be  calculated  to 

any  specified  accuracy.  This  product  basis  is  then  used  to  construct  the 

R-matrix,  or  Green’s  function,  at  the  outer  boundary  of  the  region,  R.  Matching 

to  the  asymptotic  form  of  the  wave  function  then  gives  the  S-matrix  and 

probabilities  for  rotational  transitions.  Full  details  of  the  R-matrix 

33 

calculations  will  appear  later  . 

The  laterally  averaged  potential  vQ(z)  was  modelled  by  a  Morse  potential 
and  also  by  a  variable  exponent  potential.  The  Morse  potential, 

V.M  .  -  a 

a  two  parameter  potential  (zg  does  not  affect  the  eigenvalues)  used  in  the 
previous  calculation  for  HO  on  Pt(lll),  was  used  here  primarily  for  its 
simplicity.  The  variable  exponent  potential. 

Vote)  *  +1 


(15) 


a  three  parameter  potential  Including  Eq.  (14)  as  a  special  case,  (p**>), 

is  a  flexible  form  which  has  been  used  with  considerable  success  to  fit 
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bound  states  from  atomic  and  molecular  diffractive  selective  adsorption. 

36 

An  analytic  approximation  to  its  eigenvalues  has  been  given  . 

V.  THEORETICAL  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

Calculations  were  carried  out  with  three  sets  of  potential  parameters, 

which  are  summarized  in  Table  2.  The  position  of  the  potential  minimum, 

ze,  was  chosen  such  that  the  finite  value  of  the  potential  at  z=0  is  at 

least  three  times  as  large  as  the  maximum  normal  incident  beam  energy  Ez  . 

The  three  sets  of  potential  parameters  were  determined  by  fitting  our 

diffractive  selective  adsorption  data  for  (predominantly  Jb=0)  D2  on  Ag(lll)  . 

Potential  A  is  of  the  variable  exponent  form.Eq.  (15).  The  values  of  its 

three  parameters  were  arrived  at  by  carrying  out  a  non-linear  least  squares 

fit  of  the  experimentally  determined  D2/Ag( 111)  eigenvalues  to  the  approximate 

analytical  eigenvalues.  The  parameters  for  potential  B  (also  of  variable 

exponent  form)  and  Morse  potential  C  were  determined  In  a  similar  manner, 

but  with  their  well  depths  held  constant  at  the  value  given  by  our  Debve-Waller 

20 

attenuation  measurements,  0=46  meV  .  Figure  4  shows  the  three  potentials 

A,  B  and  C  as  a  function  of  z.  Figure  5  shows  the  coupling  terms  v*(z)  as  a 

function  of  z  for  potential  A  with  8=  -0.05. 

Results  of  a  typical  scattering  calculation  are  shown  in  Figure  6.  This 

shows  the  transition  probabilities  J=0-*-J^,  J^=Q,1,2,3  as  a  function  of  the 

perpendicular  component  of  incident  energy,  E  ,  for  potential  A  with  6=  -0.05. 

zi 

The  deep  minima  just  above  threshold  in  each  channel  correspond  to  the  J+l 
resonances  seen  in  the  experimental  data.  The  theoretical  results  also  show 
weaker  maxima  and  minima  at  higher  energies,  corresponding  to  J+2  and  J+3 
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resonances.  When  a  phenomenological  dissipative  term  Is  included  in 

the  Hamiltonian,  Eq.  (11),  the  weaker  features  disappear  and  only  the 
33 

J+l  minima  remain  .  In  the  absence  of  any  dissipative  terms  in  Eq.  (11) 
the  R-matrix  calculation  preserves  unitarity,  so  that  at  energies  below  11 
meV,when  only  the  J=0  channel  is  open,  the  transition  probability  is 
necessarily  unity. 

For  each  potential  A,  B,  and  C,  a  series  of  R-matrix  calculations 

2 

were  carried  out  at  several  values  of  B.  Minimization  of  the  x  function 

V  -  £  r  v-  cf 

where  E^  are  the  experimental  resonance  energies,  their  standard  deviation 
and  Ecy  the  calculated  resonance  energies,  was  used  as  a  measure  of  the 
goodness  of  fit  of  a  given  set  of  parameters.  Sunmation  was  performed  over 
the  (HI  and  (H2  resonances  only,  since  the  calculation  does  not  provide  Ec 

p 

for  the  (HO  transition  (see  above).  Table  3  shows  the  variation  of  x  with 
3.  A  minimum  value  is  attained  for  6  between  -0.05  and  -0.10.  It  is  apparent 
that  all  three  potentials  achieve  a  comparable  fit  of  the  experimental  resonance 
data.  Slightly  different  values  of  an  optimum  3  are  obtained,  with  the  size 
of  the  estimated  value  increasing  with  well  depth.  For  potential  A,  the  best 
fit  is  obtained  with  3=  -0.05,  for  potential  B,  3  lies  between  -0.05  and 
-0.10,  and  for  potential  C,  3=  -0.10.  A  comparison  of  and  Ec^  for  potential 
A  with  3=  -0.05  is  given  in  Table  4. 
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These  calculations  were  all  carried  out  using  an  R-matrix  direct 
product  basis  of  60  translational  (sine)  functions  x  6  rotational  functions. 

The  R-matrix  boundary,  R,  was  taken  at  9.00ao.  Reference  to  Figure  4  shows 
that  the  size  of  the  residual  potential  at  the  boundary  R  is  considerably 
greater  for  the  variable  exponent  potentials  A  and  B  than  for  the  Morse 
potential  C.  This  has  a  relatively  large  effect  on  the  highest  vibrational 

i 

eigenvalues  as  shown  in  Table  5  and  also,  indirectly,  a  slight  effect  on  the 
resonance  energies  Ecv, since  the  residual  potential  should  be  negligible  at  ; '! 

the  asymptotic  matching  surface.  Calculations  for  potentials  A  and  B,with  rj  ; 

the  R-matrix  boundary  chosen  in  each  case  such  that  vq(R)  is  equal  to  the  ’ 

residual  Morse  potential  for  R=9.00ao  (-1.7  meV),  were  then  performed  for  j' 

the  best  value  of  3  from  Table  3,  3=  -0.05.  These  results  are  shown  in 

Table  6.  A  marked  improvement  in  accuracy  is  seen,  making  both  potentials  | 

A  and  B,  with  3*  -0.05,  of  comparable  accuracy  to  the  Morse  potential  C  with 

c  |  < 

its  best  value  of  3=  -0.10.  The  values  of  Ec  and  E*  for  potential  A  with  i 

V  V  I 

3=  -0.05  and  R=12.0aQ  are  listed  in  Table  4  next  to  those  previously  mentioned  j 

for  R=9.0aQ. 

VI.  DISCUSSION  j; 


Comparison  of  our  experimental  and  theoretical  resonance  energies  reveals 
that  3,  the  asymmetry  parameter,  has  a  value  which  falls  in  the  range  -0.10 
3  <  -0.05.  This  conclusion  is  independent  of  the  isotropic  potential  (A,  B,  or  C) 
which  is  used  in  the  scattering  calculation.  Further  comparison  of  these  results 
also  reveals  that  calculations  employing  either  of  the  variable  exponent  potentials 
A  or  B,  or  Morse  potential  C,  give  the  same  quality  fit  to  our  experimental 
resonance  energies.  This  relative  insensitivity  of  the  calculated  resonance 


uiUlflit 
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energies  to  slight  changes  in  the  isotropic  potential,  vQ(z),  as  compared 
with  their  sensitivity  to  the  asynwetry  parameter,  suggests  that  the 
magnitude  and  sign  of  6  can  be  well  determined  from  a  combination  of  RMSA 
and  DSA  measurements. 

The  important  and  complementary  difference  between  RMSA  and  DSA 
measurements  is  that  the  former  process  places  molecules  of  J^>0  (or  more 
specifically,  0b>vL)  into  bound  levels  of  the  physisorption  potential, 
while  the  latter  process  conserves  the  incident  J-state  (J^J^).  DSA 
measurements  can  therefore  be  used  to  obtain  the  isotropic  part  of  the 
potential  when  J^O,  due  to  the  absence  of  J-  and  m-  dependent  shifts  in 
this  case.  One  can  then  build  upon  the  resulting  vQ(z)  potential  by  making 
RMSA  measurements,  which  are  sensitive  to  the  anisotropy  of  the  potential. 

This  sensitivity  manifests  itself  as  ^-dependent  shifts  (for  J^O  there  is 
no  m-  dependence)  in  the  observed  resonance  energies.  Such  a  procedure  has 
been  followed  in  this  paper.  We  are  currently  refining  the  vq(z)  potential 

utilizing  additional  DSA  measurements  for  and  D ^  on  Ag(lll),  and  will 

20 

publish  our  final  results  elsewhere  . 

Both  the  sign  and  magnitude  of  the  asymmetry  parameter  6  are  of  interest. 
The  negative  value  of  b  found  in  these  studies  indicates  that  the  averaged 
interaction  potential  is  slightly  deeper  for  hydrogen  oriented  parallel  to  the 
surface  than  perpendicular  to  it.  Thus  hydrogen  has  a  slight  propensity  to 
reside  with  its  internuclear  axis  parallel  to  the  surface  in  the  bound  states 
of  vQ(z).  If*  a  sense  this  is  a  surprising  result.  At  large  molecule-surface 
distances  van  der  Waals  attractive  forces  dominate  the  interaction.  These 
forces  tend  to  orient  the  molecular  axis  perpendicular  to  the  surface  since 
the  polarizability  of  hydrogen  along  its  internuclear  axis,  ct1( ,  is  larger 
than  the  polarizability  perpendicular  to  it,  a^.  In  contrast  to  this,  at 


2: 


small  molecule-surface  distances  electron  repulsion  dominates  the  interaction, 

tending  to  orient  the  molecular  axis  parallel  to  the  surface.  The  two  forces, 

van  der  Waals  attraction  and  electron  repulsion,  therefore,  favor  different 

orientations  of  physisorbed  hydrogen.  This  qualitative  reasoning  is  supported 

by  Gordon-Kim  calculations  for  hydrogen  adsorbed  in  parallel  and  perpendicular 
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orientations  on  a  jellium  surface  .  The  negative  value  of  8  observed  in  our 
experiments  indicates  (at  least  in  an  averaged  sense  with  respect  to  bound 
level,  n)  that  the  orientational  anisotropy  of  hydrogen  in  the  bound  states 
sampled  here  is  more  sensitive  to  the  repulsive  part  of  the  potential. 

This  conclusion  is  contrary  to  the  assumption  made  in  most  previous 

calculations  using  potentials  having  the  form  of  Eq.  (8a),  namely  that  the 

anisotropy  parameter  8  is  determined  essentially  by  the  asymptotic  van  der  Waals 
31  32 

interaction  ’  .  However,  it  is  important  to  note  that  a  potential 

of  the  form  (8a)  with  e  treated  as  a  constant  does  not  possess  the  flexibility 
to  give  the  correct  lowest  energy  orientation  at  both  large  and  small  distances 
from  the  surface.  It  is  clear  that  the  value  of  8  obtained  here  is  really  a 
distance  averaged  quantity,  an  average  in  which  the  compensating  contributions 
from  large  and  small  distances  are  slightly  dominated  by  the  latter.  With 
regard  to  this,  it  would  be  of  value  to  perform  an  electron  energy  loss  experiment 
involving  a  low  coverage  of  physisorbed  hydrogen  on  Ag( 111)  which  would  accurately 
determine  the  m-splittings  for  the  n=Q  level,  from  which  a  level  specific  value  of  8 
could  be  derived.  It  is  our  opinion  that  treating  8  as  a  constant  will  ultimately 
prove  to  be  inadequate,  and  it  must  necessarily  be  treated  as  a  distance  dependent 
quantity  in  a  more  accurate  description  of  the  hydrogen/Ag( 111)  interaction. 

The  actual  magnitude  of  the  anisotropy  parameter,  6,  quantifies  the 
orientation  dependence  of  the  hydrogen/Ag( 111)  interaction.  As  mentioned 
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previously,  at  asymptotically  large  distances  from  the  surface  this  term 

comes  from  the  asymmetry  of  the  hydrogen  molecular  polarizability  tensor. 

31 

In  the  first  quantum  scattering  calculations  of  Wolken  it  was  assumed 

that  6  was  given  by  the  ratio  of  the  static  parallel  and  perpendicular 
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hydrogen  polarizabilities  ,  which  yields  Ba+0.24.  This  value  has  been 
used  in  mar\y  subsequent  quantum  scattering  calculations.  It  is  apprently 
not  correct.  This  has  been  recently  discussed  in  detail  by  Harris  and 
Feibelman  ,  who  took  into  account  both  the  asymmetry  of  the  molecular 
polarizability  and  the  angle  dependence  of  the  surface  response  to  a  dipole 
field,  finding  6=+0.05.  The  results  of  our  study  indicate  the  need  for  a 
still  more  comprehensive  theoretical  treatment  in  which  repulsive  interactions 
are  also  taken  into  account.  We  must  note,  however,  that  the  averaged  aniso¬ 
tropic  contributions  in  this  system  are  quite  small  as  evidenced  both  by  the 
small  magnitude  of  6  and  by  the  relatively  free  rotation  of  the  vibrational ly 
bound  HD. 

VII.  CONCLUSION 

In  this  paper  rotational ly  mediated  selective  adsorption  scattering 
resonances  have  been  used  to  determine  the  vibrational  levels  of  the  HD/Ag( 111) 
physisorption  potential  as  a  function  of  bound  rotational  state.  These  experi¬ 
mental  energy  levels  have  been  used  in  conjuction  with  exact  quantum  scattering 
calculations  to  study  the  anisotropic  part  of  the  laterally  averaged  molecule- 
surface  potential.  A  value  for  the  asymmetry  parameter  has  been  determined  which 
is  in  close  agreement  with  recent  theoretical  predictions.  Its  value  falls  in 
the  range  -0.10  -0.05. 

The  success  of  this  study  suggests  that  a  general  procedure  for  examining 
the  interaction  of  molecular  hydrogen  with  a  variety  of  substrate  materials 
has  been  developed.  Diffractive  selective  adsorption  studies  of  J=0 


or  D2  should  first  be  conducted  in  order  to  determine  the  laterally 
averaged  isotropic  part  of  the  potential.  Next,  rotationally  mediated 
selective  adsorption  measurements  should  be  carried  out  for  different 
bound  rotational  states.  Finally,  close-coupled  scattering  calculations 
can  be  used  to  analyze  the  ^-dependent  shifts  in  the  RMSA  data  in  order 
to  extract  the  orientation  dependence  of  the  molecule-surface  interaction 
potential . 

In  the  near  future  we  hope  to  refine  our  analysis  of  the  hydrogen/Ag( 111) 
potential  by  comparing  deconvoluted  rotationally  inelastic  transition 

r ; 

probabilities  with  those  calculated  using  our  R-matrix  techniques.  Additional 

diffractive  selective  adsorption  studies  are  also  planned  to  further  refine 

the  isotropic  part  of  the  potential  which  was  used  in  this  paper.  ■ 
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TABLE  1 

Experimental  Resonance  Energies  and  Bound  Levels  for  HD( J^)/Ag( 111)  (meV) 


V1 

V2 

V3 

Ez 

Ev  p^ 

Ez 

Ev  P(a) 

Ez 

Ev  P<*> 

n  3  7 

10.17 

-0.89  5 

.15 

6 

8.72 

-2.34  6 

.11 

5 

6.21 

-4.85  6 

.13 

29.67 

-3.45  7 

.23 

64.16 

-1.85  6 

.33 

4 

2.48 

-8.58  6 

.07 

25.95 

-7.17  7 

.24 

60.04 

-5.97  6 

.24 

3 

20.87 

-12.25  7 

• 

.25 

55.19 

-10.82  6 

.26 

2 

14.37 

-18.75  7 

.28 

48.56 

-17.45  6 

.27 

1 

40.21 

-25.80  6 

.30 

0 
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TABLE  2 

Parameters  for  the  three  potentials  used  to  model  the  laterally  averaged 
potential,  vQ(z). 

Potential  Form  Parameters 

D(meV)  X(aQ_1)  P  a(  ^q'1) 


A 

Variable  Exponent 

43.045 

0.6180 

4.9950 

- 

B 

Variable  Exponent 

46.000 

0.7005 

3.5312 

C 

Morse 

46.000 

. 

0.5956 

Equilibrium  Distance  zg  ■  2.3441aQ. 

Center  of  mass  displacement  for  HD,  6  *  0.233aQ. 
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TABLE  3 

2 

Variation  of  the  x  function,  Eq.  (16),  with  the  anisotropy  parameter, 
0,  for  potentials  A,  B  and  C. 


Potential 


X2(meV2) 


A 

+0.35 

842.9 

A 

+0.30 

722.6 

A 

+0.20 

427.4 

A 

+0.10 

218.5 

A 

0.00 

138.4 

A 

-0.025 

125.1 

A 

-0.05 

118.5 

A 

-0.075 

118.0 

A 

-0.10 

134.0 

A 

-0.20 

229.7 

A 

-0.30 

457.1 

A 

-0.35 

594.1 

B 

+0.20 

503.1 

B 

0.00 

141.7 

B 

-0.025 

121.2 

B 

-0.05 

110.7 

B 

-0.075 

109.4 

B 

-0.10 

111.1 

B 

-0.20 

215.8 

C 

+0.20 

597.4 

C 

0.00 

155.0 

C 

-0.025 

142.3 

C 

-0.05 

99.9 

C 

-0.075 

91.5 

C 

-0.10 

85.3 

C 

-0.15 

102.2 

C 

-0.20 

147.3 

C 

-0.30 

337.0 
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TABLE  4 

Positions  of  the  experimentally  observed  and  calculated  resonance  minima,  and 
the  corresponding  BVRR  energies  for  the  J=0+1  and  0=0+2  transition  probabilities 
shown  in  Figure  6.  Potential  A  was  used  with  8*  -0.05.  R  denotes  the  position 
of  the  R-matrix  boundary,  in  aQ.  The  bound  vibration  state  Is  labelled  by  n,  the 
rotational  state  by  J.  Experimental  resonances  are  denoted  by  Ev,  the  calculated 
resonances  by  Ecv  E'v  denotes  the  BVRR  energy,  Eq.  (12).  Values  of  Ecv  and 
E*  are  shown  for  calculations  performed  with  R=9.0ao,  and  R=12.0ao  (See  text). 


Transition 

V 

E 

V 

EC 

V 

E' 

V 

EC 

V 

E' 

V 

n,J 

R=9 .0 

R=9.0 

R-12.0 

R=12.0 

1 

14.37 

15.4 

15.46 

15.4 

15.46 

2.2 

0+1 

2 

20.87 

21.5 

21.62 

21.5 

21.62 

3,2 

3 

25.95 

26.0 

26.08 

26.0 

26.14 

4,2 

4 

29.67 

28.7 

28.87 

29.1 

29.26 

5,2 

5 

40.21 

40.3 

40.24 

40.3 

40.24 

1,3 

6 

48.56 

48.4 

48.31 

48.4 

48.31 

2,3 

0+2 

7 

55.19 

54.6 

54.49 

54.6 

54.49 

3,3 

8 

60 .04 

59.0 

58.95 

59.0 

59.01 

4,3 

TABLE  5 


Comparison  of  the  vibrational  eigenvalues  of  vq(z)  calculated  In  a  finite 
R-matrix  translational  basis  with  the  analytic  eigenvalues^  for  potentials 
A,  B  and  C.  Eigenvalues  are  shown  for  two  values  of  R,  the  R-matrlx 
boundary,  for  potentials  A  and  B.  The  vibrational  quantum  number  Is  n. 

The  translational  basis  consists  of  60  sine  functions. 

Potential  n  Analytic  eigenvalues,  en  R-matrix  eigenvalues 


R  =  9.0U 

R  =  12.00 

0 

-36.9538 

-36.949 

-36.949 

1 

-26.5775 

-26.572 

-26.572 

2 

-18.3353 

-18.346 

-18.346 

3 

-11.9879 

-12.030 

-12.030 

4 

-  7.2943 

-  7.422 

-  7.376 

5 

-  4.0109 

-  4.517 

-  4.135 

6 

-  1.8918 

-  1.212 

-  2.194 

7 

-  0.6878 

+  3.604 

-  0.590 

R  *  9.00 

R  *  12.79 

0 

-38.9176 

-38.912 

-38.912 

1 

-27.2635 

-27.261 

-27.261 

2 

-18.3965 

-18.422 

-18.422 

3 

-11.8514 

-11.926 

-11.926 

4 

-  7.1995 

-  7.397 

-  7.329 

5 

-  4.0500 

-  4.609 

-  4.226 

6 

-  2.0513 

-  1.058 

-  2.358 

7 

-  0.8925 

+  4.003 

-  1.162 

R  =  9.00 

0 

-39.8865 

-39  .  867 

1 

-28.9667 

-28.917 

2 

-19.7897 

-19.721 

3 

-12.3556 

-12.280 

4 

-  6.6644 

-  6.599 

5 

-  2.7161 

-  2.898 

6 

-  0.5106 

-  0.191 

7 

-  0.0479 

+  4.118 

a)  The  analytic  approximation  to  the  eigenvalues  of  the  variable 
exponent  potential  is  used  here  for  potentials  A  and  B. 
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TABLE  6 

2 

Values  of  x  »  Eq.  (16),  shown  for  calculations  performed  with  different 
R-matrix  values,  R,  for  potentials  A  and  B. 

Potential  6  R  X2(meV2) 


A 

-0.05 

9.0 

118.5 

12.0 

90.2 

B 

-0.05 

9.0 

110.7 

12.79 


91.5 
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FIGURE  CAPTIONS 

C13  Angular  distribution  of  the  reflected  number  density  for  HD 

scattering  from  Ag  (111)  along  the  <112>  direction.  Rotationally 
inelastic  transitions  are  labelled  Vertical  lines  denote 

the  predicted  peak  locations  for  a  monoenergetic  beam  with  =  12.7  X'1 
The  lower  section  Is  an  Ewald  representation  of  the  displayed  angular 
distribution. 

C 21  Rotationally  elastic  and  inelastic  transition  Intensities  versus 
normal  incident  energy  for  HD  Incident  along  the  <112 >  and  <101> 
azimuthal  directions.  The  vertical  lines  indicate  the  RMSA  resonance 
energies,  which  are  demonstrated  here  to  be  invariant  with  respect  to 
incident  azimuthal  orientation. 

C3]  Rotationally  elastic  and  inelastic  transition  intensities  versus 

normal  incident  energy  for  three  different  beam  energies:  40.6,  113.6 
and  153.0  meV.  The  vertical  lines  indicate  the  RMSA  resonance  energies, 
which  are  demonstrated  here  to  be  invariant  with  respect  to  incident 
momentum  parallel  to  the  crystal  surface. 

[4]  The  isotropic  component  of  the  laterally  averaged  molecule-surface 
potential,  vQ(z),  as  a  function  of  z,  the  normal  distance  to  the 
metal  surface.  Potentials  A  and  B  are  variable  exponent  potentials, 

Eq.  (15),  with  well  depths  of  43.0  meV  and  46.0  meV,  respectively. 
Potential  C  Is  a  Morse  potential,  Eq.  (14),  with  well  depth  46.0  meV. 

The  full  potential  parameters  are  given  in  Table  2. 
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[53  The  coefficients  v£(z)  for  the  Legendre  expansion  of  the  potential, 

Eq.(8b),  as  a  function  of  z,  obtained  from  potential  A  with  8-  -0.05.  j 

The  center  of  mass  displacement  from  the  HD  geometric  center  is 

5=  0.233a  .  j 
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C63  Transition  probabilities  for  rotational  excitation  of  HD  Incident  1 

on  Ag  (111),  calculated  for  potential  A,  with  6=  -0.05.  Transition  ] 

probabilities  were  calculated  at  normal  incident  energy  Increments  ;  J 

of  0.2  meV.  !  \ 
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Fig.4 


